Fractional charges and spin-charge separation in one-dimensional Wigner lattices 
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We study density response N(k, to) and one-particle spectra A(k, to) for a Wigner lattice model at 
quarter filling using exact diagonalization. We investigate these observables for models with short 
and long-range electron-electron interaction and show that truncation of the electron repulsion can 
lead to very different results. The spectra show clear signatures of charge fractionalization into pairs 
of domain walls, whose interaction can be attractive or repulsive and is controlled by the formal 
fractional charges. In striking contrast to a bound exciton in N(k,u>), we find an antibound quasi- 
particle in A(k, u), which undergoes spin-charge separation. We present a case of extreme particle- 
hole asymmetry, where photoemission shows spin-charge separation, while inverse photoemission 
exhibits an uncorrelated one-particle band. 
[ published in Phys. Rev. B 75, 125116 (2007) ] 
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I. INTRODUCTION 

In contrast to higher dimensions, where interacting 
electrons are renormalized into 'Landau quasi-particles' 
that near the Fermi surface still 'look' like electrons, 
drastic things can happen in one dimension^ The ele- 
mentary excitations are collective ones involving many 
electrons, spin and charge separate and move indepen- 
dently of each other, as seen, e.g., in photoemission 
(PES) experiments.^ Another case where many-body ef- 
fects lead to an apparent splitting of elementary particles 
is quantum number fractionalization. 3 Probably its most 
famous realization is the fractional quantum Hall effect. 4 
Likewise, Peierls distorted^ and charge-ordered^ 7 - one- 
dimensional (ID) systems with degenerate ground states 
have fractionally charged solitons or domain walls as el- 
ementary excitations. In fact, quantum number frac- 
tionalization in one and two dimensions are intimately 
related. 8 This Paper presents a study of dynamic ob- 
servables for a model showing both effects, spin-charge 
separation as well as quantum number fractionalization. 

When (long-range) Coulomb interaction is the dom- 
inant energy scale of a system, electrons try to mini- 
mize their energy by maximizing their distance and crys- 
tallize into a Wigner lattice (WL)3 Hubbard suggested 
this mechanism in the context of TCNQ charge-transfer 
salts^ where it is, however, difficult to distinguish from a 
Akp charge-density wave driven by a Fermi surface insta- 
bility. 6,7,10 As pointed out recently^ longer-range hop- 
ping changes the Fermi-surface topology in doped edge- 
sharing CuO-chain compounds like Nai+^CuO^ and 
Ca2+ y Y2_j,Cu50ior^ and this allows for a clear distinc- 
tion between the Fermi-surface independent WL and the 
charge-density wave. 

The elementary excitations of a WL consist of do- 
main walls (DWs). Their fractional charge follows from 
topological arguments^i 7 -^ and merely reflects the n-fold 
ground-state degeneracy at x = m/n filling - it is not 
related to the specific form of the interaction that gen- 
erates the charge order. After introducing the employed 



model Hamiltonian in Sec. [H] we show in Sec. IIIII that 
these fractional charges have direct physical meaning in 
a model with long-range Coulomb repulsion among elec- 
trons, because the resulting effective interaction between 
DWs is also Coulomb-like with a coupling constant de- 
termined by the fractional charges, their signs and dis- 
tance. Consistent with the formal charges, the interac- 
tion is attractive for the charge respons o 15 d 6 i 17 but re- 
pulsive for the electron addition (removal) process. We 
present and discuss charge response N(q, ui) and spectral 
densities A(k,uj) and show that the correspondence of 
formal DW charge and effective physical charge is tied to 
the long-range nature of Coulomb repulsion and absent 
for models with truncated electron-electron interaction. 

As the energy scale for spin excitations is much smaller 
than either the Coulomb repulsion or the kinetic energy, 
previous investigations considered WL formation of spin- 
less fermions. We will address the spin degree of freedom 
in Sec. IIVI and we will show how the interplay of spin 
and DW excitations leads to two different scenarios for 
the spectral density: (i) For nearest neighbor hopping, 
the coherent antibound DW excitation undergoes spin- 
charge separation in perfect analogy to an electron or 
hole added to the half-filled ID Hubbard model, (ii) For 
2nd neighbor hopping, we find striking differences be- 
tween the particle and hole channels, where an electron 
behaves like an independent particle while a hole decays 
into a spinon and a holon. 



II. MODEL AND METHODS 

We study the Hubbard- Wigner Hamiltonian 17 
H = *i I>L c m,. + h - c -) + *a J2( c L c i+2,° + h-c) 



+ J7 ^ ni^n hl + ^Vi^(ni - n)(n l+l - n) , 

i 1=1 i 
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FIG. 1: (Color online) One of the two ground states (a) with 
circles denoting holes and solid arrows spinless fermions, pair 
of domain walls (indicated by dashed vertical lines) at dis- 
tance d = 2 created by moving one electron via t\ (b), con- 
figurations with distance d = 4 between the DWs (c) . 



with nearest neighbor (NN) and next nearest neighbor 
(NNN) hopping amplitudes ti and ti. Operator c\ a 
( c i a) crea tes (destroys) an electron with spin a at site i, 

n i a = c\ a c i a and r;,, = n^f + n^j give the density, n is 
the average density. We focus on long-range Coulomb re- 
pulsion Vi — V/ ti&, but also discuss truncated interaction 
(^max = 3) to illustrate that truncation may have a strong 
impact on the excitation spectra. The NN Coulomb ma- 
trix element V is used as unit of energy. We treat chains 
up to L = 24 in the spinless case and L = 16 with spin 
using Lanczos diagonalization. 

To our knowledge, neither the single-particle spec- 
tral function A(k,oj) nor the density response N(q,u>) 
of a Wigner lattice have been studied before. We con- 
sequently treat here the most transparent case given 
by quarter filling, i.e., one electron per two sites. We 
will also analyze an effective low-energy model in terms 
of domain walls that is appropriate for small hoppings 
t\ , t% <C V, U ) 15 ' 16 i 17 and compare the results to those of 
the full model (JTJ that is formulated in terms of electrons 
or spinless fermions. 
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III. DYNAMICS OF THE SPINLESS MODEL 

At quarter-filling the ground state of (JTJ) is twofold 
degenerate for spinless fermions, and the lowest excita- 
tions are therefore given by domain walls with a formal 
charge ±e/2 i 6 ' ' i The existence of such DWs as elemen- 
tary excitations follows from topological considerations 
and only depends on the degeneracy of the ground state, 
not on the details of the Hamiltonian or the interaction 
that stabilizes the ground state. In the WL, where the 
hopping is not large enough to destroy the charge order 
(the gap is expected to vanish at t\jV ~ 0.2^2), the de- 
scription in terms of domain walls is useful, and we will 
show that the dynamics of the WL can be described by 
DWs and their interaction. DWs can be created from the 
perfectly ordered state by a NN hopping process t\, see 
Fig. Q] (a) and (b). Creating a pair of domain walls is pe- 



FIG. 2: (Color online) Density response N(q, uj) with t\ — 
—0.02V and ti — for different momenta q — nm/10, m — 
0, . . . , 10. The x show the eigen energies relative to the 
ground-state energy. 



nalized by V and costs potential energy. DWs can move 
by t\, see Fig. [1] (b) and (c), and their potential energy 
can depend on the distance d between the two DWs. In 
Refs. Ilafill the DW interaction for a model with long- 
range Coulomb repulsion has been discussed and been 
found to correspond to a Coulomb-like attraction as it 
would be expected between two physical charges of ±e/2 
each. In Ref. U3, an effective low-energy Hamiltonian in 
terms of DWs has been analyzed. 

The lowest eigen energies of (JXJ) with a small NN hop- 
ping t\ = —0.02 can be seen in Fig. [2] (x). At uj = 
and momenta k — and k — n, we see the two ground 
states. While this degeneracy at finite t\ is only perfect 
in the thermodynamic limit, the numerical data in Fig 
for a N = 20 site ring shows that the states at k = and 
k = ir are almost degenerate. The collection of states at 
lo ~ 0.5 ±0.08 have been shown to represent a continuum 
of two independent domain walls ^ and the states above 
lo ~ 0.8 mark the beginning of the four-DW continuum. 
In addition to the continua, we find an exciton corre- 
sponding to a bound 2-DW state at energies just below 
the 2-DW continuunijii 

In addition to the eigen energies, Fig. [2] shows the dy- 
namic charge structure factor 

N(q,uj) = ]T \(m\p q \fo)\ 2 5(u ~ (E m - E )) , (2) 

m 

where \m) and E m are the eigenstates and energies of the 
Hamiltonian, and |0o) is the ground state with energy 
Eq. For the perfect WL without quantum fluctuations, 
it shows only signals at uj = and momenta k — tt and 
k = 0. Spectral weight at finite- frequency is produced 
by fluctuations around perfect charge order and Fig. [2] 
reveals that most spectral weight is observed in the exci- 
ton, while the two-DW continuum at 0.4V < uj < 0.6V 
contains almost no weight. 

In edge-sharing Cu-0 chain compounds like the 
quarter-filled Na3Cu204 system, the NNN hopping t^ is 
expected to be larger than t\^. The eigenvalues and 
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FIG. 3: (Color online) Density response N(q, uj) with 
ti = —0.021^ and ti = —0.05 for different momenta q — 
nm/10, m = 0, . . . , 10. The x show the eigen energies rela- 
tive to the ground-state energy. 



charge response for t\ = — 0.02V,t2 = —0.05^ are shown 
in Fig. [3] for a 20-site ring. As explained in Ref. [13, the 
bound exciton has in this case its minimum at q = tt/2. 
At a critical value of ta the exciton gets soft and the WL 
with periodicity tt is destroyed and charge order with a 
periodicity ir/2 sets in. Again, we see that the exciton 
has large weight in N(q, u>) and should in principle be 
observable in charge response. 

Interestingly, the relevant lowest excitations of the WL 
continue to be pairs of DWs in the case of the one-particle 
spectral density 



A(k,u>) = Yl l<m + l4, T l^o)| 2 5(a; - (E+ - E )) 

m 

+ J2\(m-\c kA \4>o)\ 2 S(u> - (E-- E )) , 



(3) 



where |m + ) (\m~)) are eigenstates with eigenenergies E+ 
(E^) of the Hamiltonian with one particle added (re- 
moved). This observable is shown in Fig. 2] for a small NN 
hoping t\ = — 0.05 V and three different cases for the in- 
teraction in UJ): (i) NN repulsion only (l max = 1), where 
we see only a broad featureless continuum, (ii) slightly 
longer-range interaction with ^ max = 3, where we find a 
sharp quasiparticle below a continuum with small weight, 
and (iii) long-range Coulomb- repulsion with l maK = L/2, 
where the quasiparticle is observed above the contin- 
uum. We will now go on to show how these spectra 
- the continuum, the bound quasiparticle and the anti- 
bound quasiparticle - which have been obtained without 
any further assumption by exact diagonalization of the 
Hubbard- Wigner model, can be understood in terms of 
DWs and their interaction. 

The schematic illustration in Fig. \5\ shows that an ad- 
ditional electron again creates two DWs, this time both 
with an identical charge of — e/2. Again, their formal 
charge does not depend on details of the Hamiltonian: 
The total system has charge — e (the additional electron) 
and both DWs are equivalent, which gives each DW a 
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FIG. 4: Spectral density A(k, cj) for spinless fermions for (a) 
truncated Vi with l m ax = 1, (b) truncated V; with l ma x = 13 
and (c) long-range Vi with ti = —0.05V, ti = 0. Insets show 
the density of states. 
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FIG. 5: (Color online) One of the two ground states (a) with 
circles denoting holes and solid arrows spinless fermions, pair 
of domain walls (indicated by dashed vertical lines) at dis- 
tance d = 1 created by adding one electron (b), configurations 
with distance d — 3 between the DWs mediated by ti(c). 
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FIG. 6: (Color online) Potential Ed between two domain walls 
at distance d, see Fig. [5] for long-range electron-electron in- 
teraction (l max = L/2, V), effective DW interaction Q (*), 
and for truncated electron-electron interaction (7 m ax = 1, x 
and I m ax = 3, o). 



formal charge — e/2. The interaction Ed between two do- 
main walls at distance d can be obtained by calculating 
the potential energy for a configuration where they are d 
sites apart (as shown in Fig. [5] for distances d = 1,3) by 
setting all hoppings to zero. The resulting interaction is 
depicted in Fig. [6] both for long-range and for truncated 
electron-electron interaction. If the Hamiltonian fTJ) has 
only NN Coulomb repulsion Z max = 1, the energy (x in 
Fig. O does not depend on d, and the DWs are therefore 
independent. For Z max = 3, the DWs are independent 
for large distance and attractive for small d. This means 
that their interaction has the opposite sign from the one 
expected from their formal charges of —e/2 each. For the 
long-range case / max = L/2 (V in Fig. [6]), however, the 
asymptotic behavior of the interaction between two DWs 
decays like their inverse distance d: 



E d /V~l/2 + e 2 /(4d) 



(4) 



Hence, the interaction between the DWs is Coulomb-like 
and the prefactor is given by the two fractional charges 
of —e/2. The asymptotic relation (* in Fig. O gives an 
excellent approximation already for rather small DW dis- 
tances d > 1. That (formal) fractional charges are the 
elementary excitations in charge ordered chains can be 
derived from the ground-state degeneracy alone^ with- 
out any reference to the particular form of the interac- 
tion stabilizing the degenerate ground states. In the case 
of long-range electron-electron interaction, and only in 
this case, the formal fractional charges have a very direct 
physical meaning: Their interaction corresponds exactly 
to that expected for 'half-electrons'. 

For the description of the photoemission process of a 
spinlcss fcrmion added to the perfect WL an effective 
low-energy Hamiltonian in terms of DWs can be obtained 
that contains hopping t\ and the potential energy Ed- A 
reader interested in details for the following short deriva- 
tion should consult Ref. [l7|, where an analogous treat- 
ment is given for the case depicted in Fig. [TJ Starting 



from the ground state with p = 
1 



V2 



(|o t o t o t •• •) + 1 T o T o T o.. .)) (5) 



and adding an electron with momentum fc, we arrive at a 
state \Tp d= \,k) with momentum k where the DW centers 
have distance d = 1 (see Fig. [SJb)) 
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o TTT o T 
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(6) 



where the — signs are due to a Fermi sign obtained by 
moving the creator through an existing electron. Once 
the two DWs have been created at distance d = 1, they 
can move via NN electron hopping t±: 

L 

Ti\tpd=i,k) = h X!( c r,t c »'+i,T + n - c -) \^d=l,k) = 



= (e ik (\ t o TT o T • • ■ ) + I TT o T ■ • • TT o>) 



2L 

-e 2lkr (\ TT o TT o . . . ) 
-e 3ife (| T o TT o T •■■) 

+...) = 

2ti sink\ip d =3,k) , 



|o TT o T o. 
|o TT o TT • 



■ TT)) 
•)) 

(7) 



where \ipd=3,k) denotes the state with momentum k and 
distance d = 3 between the DWs (a complex phase has 
been absorbed into its definition). Repeating step |7]), 
d can grow further, which each t\ process increasing (or 
reducing) d in steps of two, see Fig. [5] (Additionally, 
t\ can introduce more DWs, as in Fig. [T]Jb), but such 
processes cost energy and will be neglected in this low- 
energy analysis.) 

We now arrive at the effective DW Hamiltonian 7Ydw 
in terms of d and k. 



T~Cdw = 



/ E x fi(fc) ... \ 

ti(fc) E 3 h(k) ... 

fi(jfe) E 5 tt(k) .. 

V ; J 



(8) 



where the diagonal contains the potential energy E d 
and momentum k enters the effective hopping ^i(fc) = 
2t\ sin k. 

We are interested in the spectral density (O, which can 
be written as 



A(k,u) 



1 



■~3^ |c feT - 
tt cj — Hbw + i0 



(H-E o )+i0+ fe < T ' 



(9) 
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FIG. 7: Spectral density obtained from the effective two- 
DW Hamiltonian © for ti = -0.05V, t 2 = 0. (a) For 



non-interacting DWs with Ed = const, see Fig. 4(a) for 
the spectral density for the corresponding case i max = 1 in 
the electron Hamiltonian (JlJ. (b) Attractive DW interac- 
tion with Ei = 0.83V, E d> - i — V, resulting from i max = 3 
in <[TJ) , compare to Fig. |4(b)| (c) Repulsive DW interaction 
with Ei — 0.7V, Ed>3 = 0.5V as an approximation to the 
DW potential of the long-range electron Hamiltonian with 
^max — L/2, the pertaining spectral density is seen in Fig. |4(c)| 
The 6 peaks of bound states have been slightly broadened 
and their height has been cut for better visibility of the 2 
DW-continua. 



i.e., we need the (1, 1) element of the inverse of the matrix 
lo — 7Yd\Vj which can be written as a continued fraction 



Hdw)i,i = 



lo 



ti(fc) 2 



(10) 



This fraction can be easily evaluated for the distance in- 
dependent DW potential Ed — V that is found, if only 
NN Coulomb repulsion is kept in the electronic Hamilto- 



nian (JTj), see Fig.[6l We then obtain 



A\ — (lo — 7xdw)i i 



E 



1 



where E is denotes the distance-independent DW po- 
tential, here E = V. The resulting spectral density 
has only a regular part and no singularities (except for 
k = 0, 7T, where the hopping t\(k) vanishes), and is shown 
in Fig. 7(a) The incoherent two-DW continuum at ener- 
gies V — 4ti sin k < lo < V + 4ti sin k results from the in- 
dependent (since Ed = E = const) movement of the two 
DWs, and c orresponds to the spectral density depicted 
in Fig. 



4(a) which was calculated for the full fermionic 



model yj) with only NN repulsion (i„ 
Next, we move to the case Z ma x 



1). 



3, which leads to 
a short range attraction between the DWs, see Fig. [6] 
Since the potential is constant for d > 3, we can use the 
result obtained above for Ed = const in order to arrive 
at 



A, = 



1 



Ei-h(k)*Ai ' 



(12) 



which is shown in Fig. |7(b)| In addition to the contin- 
uum, we now find a pole where the denominator of (|12| 
vanishes. It results from the DW attraction and corre- 
sponds to a bound quasiparticle with dispersion 



e(k) = E 1 



h(k) [ 



E\ — E 



(2iisinfc) 2 
Ei - E 



(13) 



where E = V is the potential at distances d > 3 and 
Ei ~ 0.83V the potential at d = 1. For long-range 
electron-electron repulsion, we simplify the resulting DW 
potential to E\ — 0.7V, E = 0.5V and can then apply 
(|12[) and (fT3|) . However, the quasiparticle is now an anti- 



bound state above the incoherent continuum, see Fig 7(c) 



This highly unusual situation is also observed in the spec- 
tra of the full model ( see Fig. 4(c) ) obtained by exact 
diagonalization. As we have seen here, it results from 
the DW repulsion and is thus a direct consequence of the 
long-range nature of the Coulomb repulsion absent from 
short-range models. 

The only feature of the spectra for the full model (Q]) 
that is not reproduced in the simplified DW analysis is 
the transfer of spectral weight seen in Fig. [4] with more 
weight in photoemission (PES) at k = and more in 
inverse photoemission (IPES) at k = n. This transfer 
is caused by fluctuations around perfect charge order in 
the ground state induced by ti, see Fig. [T] and will be 
discussed elsewhere. 

Motivated by experimental data indicating that NNN 
hopping t-x is larger than ti in Na3Cu204jii we now 
turn to the spectral density for t% ^ 0. Via NNN hop- 
ping, an electron inserted into an empty site of the WL 
can hop over the occupied sites and move along the 
empty ones like a free fermion, and the same holds for 
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FIG. 8: Spectral density A(k,u>) for ti = -0.02V, t 2 = 
—0.05V. (a) For the spinless-fermion Hamiltonian ([TJ with 
long-range interaction Z max = L/2 obtained by Lanczos di- 
agonalization via ©. (b) For the effective DW-Hamiltonian 
© obtained by dHJ) with E x = 0.7V + 2t 2 cos2fc, E = 0.5V. 
Insets show the density of states. 



a hole inserted into an occupied site. For t\ = 0, we 
therefore obtain an independent-electron-like dispersion 
with e(fc) = Vln2 + 2i 2 cos2fc (not shown). The dis- 
persion is the same in PES and IPES, which results in 
an indirect band gap with the highest occupied state at 
k = 7r/2 and the lowest unoccupied ones at k = 0, ir. 
Finite NN hopping t\ ^ can be taken into account 
just like in the t 2 = case discussed above by setting 
Ei = Vln2 + 2i 2 cos2fc in The resulting spe ctral 

density of the effective DW model is shown in Fig. |8(b)| 
and the corresponding results for the fermion Hamilto- 
~ Both for PES and for IPES, we see 



nian (p} in Fig. 8 (a 



the 2-DW continuum as well as a quasiparticle dispersion 
~ cos 2k. But we observe a certain particle-hole asymme- 
try: For IPES, the continuum has small weight and the 
width of the dispersion is 4i 2 as for a free electron. For 
PES, the continuum is stronger and almost mixes with 
the quasiparticle at momenta k ~ 7r/2, which somewhat 
reduces the band width for the quasiparticle. We will see 
in the next section that the particle-hole asymmetry is 
strongly enhanced for electrons with spin. 



FIG. 9: (Color online) An electron added to the perfect WL, 
as Fig. [5] but for electrons with spin. In addition to the 
two DWs, an electron now creates a spinon (denoted by the 
unfilled circle/ellipse). As in the spinless case, the DWs can 
move via ti, and the spinon can move via spin-flip processes. 




FIG. 10: (Color online) Spectral density A{k,u)) with spin 
for t x = -0.05V, t 2 = 0, U = 4V, L = 16. The right-hand 
inset shows a blow-up of IPES with a higher energy-resolution, 
the dashed and solid lines indicate spinon and 'anti-holon' 
branches (guides to the eye). The inset on the left shows the 
density of states. 



IV. DYNAMICS OF ELECTRONS WITH SPIN 

Having understood the dynamics of the spinless sys- 
tem, i.e., of the charge sector, we now include the spin. 
An inserted electron or hole then produces a spinon in 
addition to the two DWs, see the cartoon in Fig. [5] il- 
lustrating the possible processes in the WL with spin. 
Again, the dynamics of the charge sector are determined 
by DW hopping via t\. The spinon can move by spin- flip 
processes ~ J 2 = t^/U or ~ J\ = t^/ll, a magnetic scale 
much smaller than the hopping scale. In this section, 
we will discuss how the spectral density indicates spin- 
charge separation similar as in the half-filled ID Hubbard 
model. 

At first, we will focus on NN hopping and compare 
the spectral density for t-y — —0.05V, t% — with spin 



(Fig. [TO]) to the one without spin (Fig. 4(c) I. The 2-DW 
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FIG. 11: Comparison of the spectral density for Wigner and 
Hubbard models, (a) Blow-up with higher energy resolution 
of IPES for the Wigner model (L — 16 sites) as in Fig. [10] 
The dashed and solid lines are guides to the eye indicating 
spinon and 'anti-holon' branches, (b) IPES for the half-filled 
ID Hubbard model (L = 8) with NN topping t = 1 and 
U = 200t. 



While this may not be so obvious on the present short 
8-site ring, the anti-holon with width ~ t and the spinon 
~ J = 4t 2 jU <C t are clearly visible in the thermody- 
namic limit. 20 The strong similarities seen in the spectral 
densities Fig. I 111 for the Wigner and the Hubbard models 
let us conclude that the Wigner model likewise shows a 
convolution of spinon and anti-holon bands; the role of 
the anti-holon is now taken by the antibound 2-DW state 
that is the elementary collective excitation of the charge 
sector, see Fig. [TST c). Since the antibound 2-DW state 
already has a periodicity of 7T, see (|13|) . the momentum 
range k — to k = it/2 corresponds to k = to k = it 
in the Hubbard model with NN hopping, where the one- 
particle dispersion has periodicity 2tt. The width of the 
convolution is determined by that of the antibound 2-DW 

state, hence ~ ^^y-, see (|13p ; the spinon has almost 
no dispersion, because ~ J\ <C t±. The small differences 



between the Figs. 11(a) and 11(b) are due to weak inter- 
actions with the 2-DW continuum in the former case. 



(a) 






DWBS \ 



(c) 



FIG. 12: (a) Charge-spin separation in ID Hubbard model: 
Electron (solid line) decays into holon and spinon; (b) Wigner 
lattice at quarter-filling: Decay of electron into 2 DWs and 
a spinon; (c) Long-range Coulomb-repulsion of DWs gener- 
ates a DW (anti) bound state (DWBS). This scattering state 
behaves like a holon. 



continuum is broader with spin, e.g. at k = and/c = tt, 
where the width of the continuum shrinks to zero in the 
spinlcss model. A more fundamental change is apparent 
in the antibound quasi-particle of the spinlcss case: It 
has evolved into a narrow structure with high spectral 
weight comprised of several peaks. Figure 11(a) shows a 



blow-up with higher energy resolution of this structure 
(in IPES) for k = to k — tt/2. When we compare 
this blow-up to IPES for the usual half-filled ID H ub- 
bard model (t = 1, U = 200i) shown in Fig. |ll(b)| we 
observe remarkably similar structures. It is well known 
that an electron in the Hubbard model decays into a 
spinon and an anti-holon, see Fig. [T27 a). and that IPES is 
given by a convolution of spinon and anti-holon branches. 
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FIG. 13: (Color online) Spectral density A{k,u)) with spin 
for (a) ti = 0,* 2 = -0.051/, (b) ti = -0.02V, t 2 = -0.05V, 
U — 41/, L = 16. Insets on the left show the density of states. 
In (a), spinon and holon branches in PES are indicated by 
dashed and solid lines (guides to the eye), the dash-dotted 
fine in IPES gives the one-particle dispersion E(k) = V In 2 + 
2t2 cos 2k. 



We now turn to the spectral density for t% ^ 0, at first 
choosing t\ = for simplicity. The ground state is then 
given by the perfectly charge ordered WL, an electron 
added in IPES goes into an empty site and can move 
freely on the empty sublattice, without any interaction 
with the occupied sublattice (as t\ =0). Consequently, 
IPES shows a one-particle-like tight-binding band with 



1 



8 



dispersion E(k) — Fin 2 



2^2 cos 2k just like in the spin- 
The situation for 
however, fundamentally different: The hole 



less model, see u> > in Fig. 13(a) 
a hole is 



goes into the occupied sublattice, which is a half-filled 
Hubbard chain, and it therefore separates into spin and 
charge. The resulting spinon and holon branches in PES 
can be seen for lu < in Fig. 13(a) (Again, the spectrum 



was found to agree to the one for the Hubbard model 
on eight sites. In this case even without any deviation, 
because there is no 2-DW continuum.) In this case of ex- 
treme particle-hole asymmetry, one and the same observ- 
able, the spectral density, shows both pure one-particle 
behavior (in the particle sector) and strongly correlated 
behavior (in the hole sector) . If both and t\ are active, 
the inserted particle (hole) can interact with the occupied 
(empty) sub-lattice, see the spectral density in Fig. |13(b)| 
t\ -hopping processes now induce an additional incoher- 
ent 2-DW continuum, and we aquire incoherent weight in 
IPES. However, the strong particle-hole asymmetry per- 
sists: As in the spinless case shown in Fig. [5J we see that 
the quasiparticle in IPES remains largely intact and that 
- surprisingly - the peaks furthest from the Fermi energy 
remain sharpest. 



V. SUMMARY AND CONCLUSIONS 

To conclude, we have investigated the dynamics of a 
quarter-filled Hubbard- Wigner model. Apart from be- 
ing the most transparent situation on a Wigner lattice, 
quarter-filling is also appropriate for the edge-sharing 
chain compound Na3Cu2 04. We find that electrons de- 
cay into a spinon and two domain walls as sketched in 
Fig- HUM. 



In the absence of the spinon, i.e., for the spinless model, 
we have investigated the effects of a truncated Coulomb 
interaction vs. fully long-range interaction in Sec. IIIII 
By comparing dynamic observables of the Wigner model 
(Q]) to those of an effective DW Hamiltonian, we have 
shown that the DWs and their interaction clearly man- 
ifest themselves in observables and are thus accessible 
to experiment. In the case of truly long-range electron- 
electron interaction, the formal fractional charges can be 
given a direct physical meaning that is consistent for both 
one-particle spectra A(k,uj) and two-particle dynamics 
N(q,uj). In models with truncated interaction, we still 
find formal fractional charges, but their interaction does 
no longer directly correspond to their formal charge. This 
difference has a strong impact on A(kuj): The DW repul- 
sion due to long-range electron-electron interaction leads 
to a highly unusual anti-bound quasiparticle vs. a more 
conventional bound quasi-particle observed for truncated 
electron-electron interaction. 

We have finally analyzed the Hubbard- Wigner model 
for electrons with spin, where we find the signatures of 
both charge fractionalization and spin-charge separation. 
Despite its composite nature, the antibound quasiparti- 
cle undergoes spin-charge separation reminiscent of an 
electron in the ID Hubbard model, see Fig. [T2Tc). Ex- 
perimentally, £2 > t\ could be more relevant^ and this 
case shows striking particle-hole asymmetry: For t\ = 0, 
an added electron behaves like an independent particle 
while a hole shows strongly correlated behavior and de- 
composes into spinon and holon. Even for finite t\, this 
behavior persists to some extend. 

We thank A. M. Oles, J. Unterhinninghofen and 
O. Gunnarsson for careful reading of the manuscript and 
for helpful suggestions. 
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